Superfluidity of "dirty" indirect excitons in coupled quantum wells 
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The theory of what happens to a superfluid in a random field, known as the "dirty boson" problem, 
directly relates to a real experimental system presently under study by several groups, namely 
excitons in coupled semiconductor quantum wells. We consider the case of bosons in two dimensions 
in a random field, when the random field can be large compared to the repulsive exciton-exciton 
interaction energy, but is small compared to the exciton binding energy. The interaction between 
excitons is taken into account in the ladder approximation. The coherent potential approximation 
allows us to derive the exciton Green's function for a wide range of the random field strength, and 
in the weak-scattering limit CPA results in the second-order Born approximation. For quasi-two- 
dimensional excitonic systems, the density of the superfluid component and the Kosterlitz-Thouless 
temperature of the superfluid phase transition are obtained, and are found to decrease as the random 
field increases. 
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Superfluidity in a system of spatially indirect ex- 
citons in coupled quantum wells (CQW) has been 
predicted by Lozovik and Yudson^^ and several sub- 
sequent studiesSiSi^iSiSiLSiSiiSiiiiiS have predicted that 
this should be manifested as persistent electric cur- 
rents, quasi- Josephson phenomena and unusual proper- 
ties in strong magnetic fields. In the past ten years, 
a number of experimental studies have focused on 



this goal 
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The coupled quantum well 



system is conceptually simple: negative electrons are 
trapped in a two-dimensional plane, while an equal num- 
ber of positive holes is trapped in a parallel plane a 
distance D away. One of the appeals of this system is 
that the electron and hole wavefunctions have very lit- 
tle overlap, so that the excitons can have very long life- 
time (> 100 ns), and therefore they can be treated as 
metastable particles to which quasiequilibrium statistics 
apply. Also, when D is large enough, the interactions be- 
tween the excitons are entirely repulsive, so that the ex- 
citon gas is stable against collapse (see, e.g., Ref.L21 22]). 
In general, the theory of these systems has neglected 
the role of disorder. In real experiments, however, disor- 
der plays a major role, arising from impurities, boundary 
irregularities, and fluctuations of the alloy concentration 
of the epitaxial layers. Although the quality of quantum 
well structures has improved dramatically over the past 
decade, e.g., the inhomogeneous broadening linewidth 
of typical GaAs-based samples has been improved from 
around 20 meV to less than 1 meVji^ the disorder energy 
is still not negligible compared to other energies. At a 
typical exciton density of 10^° cm~^, the interaction en- 
ergy of the excitons is approximately Airh^an/MD ~ 1 
meV (a being the characteristic radius of a 2D exciton, M 



the total mass, and n the number density of excitons) iS£ 
Typical thermal energies at liquid helium temperatures 
are ksT = 0.2 — 2 meV. On the other hand, the typical 
disorder energy of 1 meV is low compared to the typical 
exciton binding energy of 5 meV. 

Earlier studies of disorder in exciton systems include 
theory of the transport properties of direct and indi- 
rect excitons and magnetoexcitons in random fields, S^ 
the influence of various random fields on excitonic and 
magnetoexcitonic absorption of light )S^ and Anderson 
localization of excitonsiSi In two-dimensional systems, 
the excitonic interaction in the Bogoliubov approxima- 
tion is valid only at non-physically low densities be- 
cause of the divergence of the two-dimensional scattering 
amplitude in the Born approximationi2^ Therefore, the 
ladder approximation must be used at low densities to 
treat properly the interaction between two-dimensional 
excitonsi2S*2& In this letter we study the case of a ran- 
dom field which is not necessarily small compared to the 
dipole-dipole repulsion between excitons. The coherent 
potential approximation^"* (CPA) allows us to derive the 
2D indirect exciton Green's function for a wide range of 
the random field, and in the weak-scattering limit CPA 
results in the second-order Born approximation. 

As mentioned above, fluctuations of the thickness of 
the quantum wells, which arise during the fabrication 
process, impurities in the system, and disorder in the al- 
loy of the barriers can all lead to the appearance of a 
random fleld. Of these, spectral analysis of the exciton 
luminescence shows that alloy disorder, with a character- 
istic length scale short compared to the excitonic Bohr 
radius of around 100 A, plays the most important role. 
In our model, the random potential ^(r) acting on an 



electron and hole is considered to be a Gaussian noise, 
such that 



{V{r)V{r'))^gd{r-r'), 



{V{r)) = 0. (1) 



An electron is subjected to the potential Ve = ae^(r) 
and a hole to Vh = —ahV{r), where a^ and ah are con- 
stants. We consider the characteristic length of the ran- 
dom field potential I to be much smaller than the aver- 
age distance between excitons r^ ~ l/y^Tm (I <^ l/y^Tm, 
where n is the total exciton density). Therefore, in or- 
der to obtain the Green's function of the excitons with 
dipole-dipole repulsion in the random field, we obtain the 
Green's function of a single exciton in the random field 
(not interacting with other excitons), and then apply per- 
turbation theory with respect to dipole-dipole repulsion 
between excitons, using the gas of non-interacting exci- 
tons as a reference system. 

The momentum-frequency domain Green's function of 
the center of mass of the isolated exciton in the random 
field at T = within the coherent potential approxima- 
tion (CPA) is given by— (here and below h = I) 



G<"'(P, 



1 



Lu - ea{p) + ^1 + iQ{p, t^) ' 



(2) 



where /i is the chemical potential of the system, and 
Snip) — p^ /2M is the spectrum of the center mass of 
the exciton in the "clean" system (M = r^e -I- m^ is the 
mass of the exciton; rUe and mu are the electron and hole 
masses, correspondingly). The function (5(p,w) is deter- 
mined by the effective random field acting on the center of 
mass of the exciton. For zero random field, Q(p, uj) -^ 0. 
li g -^ El, {Ef, is the binding energy of an exciton) , the 
function Q{p, uj) in the second-order Born approximation 
can be expressed through the ground-state wavefunction 
of an indirect exciton. Analytical results for (5(p, a;) may 
be derived for D ^ a and me — rrih or to/,, ^ rrie and 

ah ( ^ ) <C ag. We present in this Letter the simplest 

case rrie — mh- The result, for the case D ^ a and 
TOe = mh, is 



Q(p,w) 
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Here p = (Sa)!/''!?^/'', with a = e/Ame^he^ (e is the 
dielectric constant of the material, rue-h is the reduced 
mass of exciton, and e is the electron charge), and Jq{z) 
is a Bessel function of the first kind. 

At small densities n [np^ <^ 1), the system of indi- 
rect excitons at low temperatures is a two-dimensional 
weakly nonideal Bose gas of dipoles with dipole moments 
d normal to wells (d ~ eD). The distinction between 
excitons and bosons manifes ts itse lf in exchange effects 
(see, e.g., Refs.0i3, and ^2S\). These effects are 
suppressed for excitons with spatially separated e and 



/i in a dilute system (np^ <C 1) at large D {D ^ a), 
because at large D, the exchange interaction in the spa- 
tially separated system is suppressed relative to the e—h 
system in a single well due to the smallness of the tun- 
neling exponent T ^ exp[— (D/2a)^'^] originating from 
the dipole-dipole interaction. Hence, when D ^ a, ex- 
change phenomena, connected with the distinction be- 
tween excitons and bosons, can be neglected, and there- 
fore, the system of indirect excitons in CQWs can be 
treated by diagrammatic techniques developed for boson 
systemsSS,. Two indirect excitons in a dilute system in- 
teract as U{R) = e^D^/{eR^), where R is the distance 
between exciton dipoles along quantum well planes. 

At the characteristic frequencies and momenta which 
give the greatest contribution to the Green's function in 
the ladder approximation, the function Q{p,uj) can be 
approximated by the constant (see Eq. ^) Q{p — 0,uj — 



0) = Q = (ae 



ahj 



^Mg/(167r*). As a resuh, the first 



order approximation in the amplitude of scattering of the 
isolated pair of bosons /o the two-particle vertex F can 
be reduced to /o pertinent to a clean syste m^^i^^ , i.e. F 
does not depend on Q in the first approximation. As a 
result, the chemical potential /i has the same form as for 
the pure systemS^: 
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We have the condensate Green's function D{p,iojk) 
n(o)/ ■ N -^(27r)^no(5(p) 

where Uq is the density of Bose condensate. Since at 
small temperatures (n — nQ)/n <^ 1, according to the 
ladder approximation^^ we use n instead of hq below. 
G{p,iLUk) and F(jp,iujk) are the normal and anomalous 
Green functions of the overcondensate: 

^, . ^ iujk + £o{p) + P + -iQ 
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where £o{p) is the spectrum of noninteracting excitons; 
for small momenta p <ti p the excitation spectrum of the 
interacting excitons, e{p), is acoustic: e{p) = c^p, where 



Cs = y y/p'^ — Q'^/M is the velocity of sound. 

We calculate the density of the normal component 
nn{T), using the Kubo formula^^ and the total single- 
particle Matsubara Green's function of the indirect exci- 
tons. We obtain 
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Here N is the total number of particles, and n° is the 
density of the normal component in a pure system with 
no impurities: 

^ 3_ [ ^P 2 ^"o(p) 

""~ 2M J (27r)2^ de ' 



(8) 



where no{p) = (e.^'^py^ — 1)~^ is the distribution of an 
ideal Bose gas of temperature excitations. 

The first term in Eq. (0, which does not depend on Q, 
is the contribution to the normal component due to scat- 
tering of quasiparticlcs with an acoustic spectrum in an 
ordered system at T 7^ 0. In a two-dimensional system, 



3C(3) T^ 
1-K di(n,Q)M' 



(9) 



where C,{z) is the Riemann zeta function (C(3) ~ 1.202). 
The second term in Eq. iQ is the contribution to the 
normal component due to the interaction of the particles 
(excitons) with the random field, which, when explicitly 
evaluated, yields: 



nQ 



2Mcl{n,Q\ 



(10) 



The density of the superfluid component is n^ = n ~ 
Tin. From Eqs. and H10|) we can see that increasing 
the random field decreases the density of the superfluid 
component. 

In a 2D system, superfluidity appears below 
the Kosterlitz-Thouless transition temperature T^ = 
7rns/(2M)^, where only coupled vortices are present. 
Using the expressions Eqs. © and (fTn|l for the density Ug 
of the superfluid component, we obtain an equation for 
the Kosterlitz-Thouless transition temperature T^. Its 
solution is 



Tr. = 
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Here T^ is an auxiliary quantity, equal to the temperature 
at which the superfluid density vanishes in the mean-field 
approximation ng(T^) — 0: 
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(12) 



furthermore, in Eqs. Hll() and (|12|l . 



n = n — 



nQ 

2Mc1 ' 



(13) 



The dependence of the Kosterlitz-Thouless transition 
temperature T^ as a function of the total exciton density 
n for different Q, obtained from Eq. (|11|) . is presented in 
Fig. 1. It can be seen in Fig. 1 that the random field de- 
creases the Kosterlitz-Thouless transition temperature. 
Fig. 2 shows the dependence of the KTS transition tem- 
perature on the random field parameter Q. 

This work shows that although the random fleld de- 
pletes the condensate, Kosterlitz-Thouless superfluidity 
should still be possible in a system of spatially indirect ex- 
citons. For realistic experimental parameters taken from 
photoluminescence line broadening measurementsi^, the 
disorder energy is approximately 1 meV, which implies 
Q = 0.05 in the units of the plots here, which implies a 
KTS transition temperature of T = 12 K at a density 
of 3 X 10^" cm^^. This is still well below the density 
n = l/p{D)'^ at which breakdown of the boson picture 
of excitons can occur, which for a typical electron-hole 
separation of 120 A , i.e. Z? = 3 in the units used here, 
is approximately 4 x 10^^ cm~^. On the other hand, for 
larger values of Q, which corresponds to sample quality 
of just a few years ago, the critical density can become 
quite close to this limit. 
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Captures to Figures (1-2) 

Fig.l. Dependence of temperature of Kosterlitz- 



Thouless transition Tc = Tc{n) (in units of Ry* 



Vea; 



a — (4TOe-/ie^/e)~^) on the exciton density n (in units of 
a~^) at the interwell distance Z) = 3 (in units of a); at 
different random fields Q (in units of Ry): Q = - solid 
curve; Q = 0.05 - dotted curve; Q = 0.1 - dashed curve; 
Q = 0.2 - dashed-dotted curve. 

Fig. 2. Dependence of temperature of Kosterlitz- 
Thouless transition Tc = Tc{Q) (in units of Ry* = e^/ea; 
a = (4TOe-he^/e)^^) on the random field Q (in units of 
Ry) at the interwell distance Z) = 3 (in units of -Ry); at 
the different exciton densities n: n = 0.005 - solid curve; 
n = 0.007 - dotted curve; n = 0.01 - dashed curve. 
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Figure 2 
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